Phase diagram of 4 He adsorbed on graphite 
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We present results of a theoretical study of 4 He films adsorbed on graphite, based on the contin- 
uous space worm algorithm. In the first layer, we find a domain-wall phase and a (7/16) registered 
structure between the commensurate (1/3) and the incommensurate solid phases. For the second 
layer, we find only superfluid and incommensurate solid phases. The commensurate phase found in 
previous simulation work is only observed if first layer particles are kept fixed; it disappears upon 
explicitly including their zero-point fluctuations. No evidence of any "supersolid" phase is found. 
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I. INTRODUCTION 

Adsorption of helium on a graphite substrate is still 
the subject of many experimental and theoretical studies; 
although this subject is almost four decades old, it has 
lately been enjoying a resurgence of interest in connection 
with the study of a possible supersolid phase of matter. 

Due to the strong attraction to graphite, helium forms 
up to seven distinct layers above the substrate;^ each 
layer being a realization of a quasi-two dimensional sys- 
tem. Several types of phases result from the interplay 
between the interaction among helium atoms, and their 
interaction with the substrate, including fluid, commen- 
surate and incommensurate solid phases. Crowell and 
Reppy raised the possibility of a supersolid phase in the 
second layer, from the anomalous behavior of the pe- 
riod shift in torsional oscillator experiments. 3 Supersolids 
exhibit simultaneously crystalline order and frictionless 
flow in a single homogeneous phase, and have attracted 
increasing interest since the observation of nonclassical 
moment of inertia in solid He by Kim and Chani^ It 
has been proposed, based on a number of fundamental 
arguments, that a commensurate perfect single crystal 
of 4 He ought not be supersolid, 5 a prediction supported 
by a number of computer simulations^' - However, a su- 
persolid phase exists for bosonic models on a triangular 
lattice^ and one may thus speculate about the possibility 
of supersolids on substrates^ 

For the first adsorbed helium layer, there exist pre- 
ferred adsorption sites, located above the centers of the 
hexagons formed by the carbon atoms on the graphite 
surface. A commensurate phase at filling 1/3, as well 
as an incommensurate solid phase are clearly observed 
in neutron diffraction experiments ; 10 ' 11 ! 12 heat capac- 
ity measurements ; 13 ! 14 ' 15 and in numerical simulations. 16 
The identification of the phases occurring between the 
two solids is still uncertain^ Several types of domain- 
wall phases have been predictedji^ but none unambigu- 
ously observed. 

The second layer is known to exhibit a gas, a su- 
perfluid and an incommensurate solid phase, as shown 
by heat capacity measurement o 15 ' 19 and neutron diffrac- 



these two phases, Greywall and Busch conjectured a com- 
mensurate solid with a V7 x \fl partial registry with 
respect to the first layer, based on their heat capacity 
measurements jiii 5 - At the same filling, a commensurate 
solid phase was observed in path integral Monte Carlo 
(PIMC) simulations^ - in that study, however, for com- 
putational convenience first layer particles were treated 
as classical, i.e., held fixed in space at their T=0 equi- 
librium position. Because of the relative weakness of the 
adsorption potential experienced by second layer atoms, 
it is plausible that the explicit inclusion of zero-point mo- 
tion of first layer particles, possibly leading to a further 
weakening of the attraction, may qualitatively alter the 
picture. Moreover, no prediction has yet been made the- 
oretically regarding the existence of a possible supersolid 
phase in the second layer. 

In this paper, we study the low temperature phase dia- 
gram of the first and second layers of 4 He on graphite, by 
means of state-of-the-art computer simulations in which 
quantum zero-point motion of helium atoms in both the 
first and second layers is fully included. 

Our computed phase diagram is summarized in Fig. [TJ 
In the first layer we find a striped phase, a hexago- 
nal domain-wall phase, and a hexagonal commensurate 
structure (at filling 7/16) between the two crystals. In 
the second layer, we do not observe a commensurate solid 
phased in contrast to the predictions by Greywall and 
Busch. 1,15 Rather, as coverage is increased the system 
goes through coexistence of liquid and gas phases, a (su- 
perfluid) liquid and an incommensurate crystal. 

We show that the commensurate phase observed in 
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FIG. 1: Schematic phase diagram of the first (a) and second 
(b) layers, as a function of two-dimensional layer density. 
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Ref. is merely a consequence of a the neglect of zero- 
point motion of first layer 4 He atoms, and that such a 
phase disappears if this approximation is removed, i.e., if 
zero-point motion of first layer particles is included in the 
simulation. In no case is a supersolid phase of 4 He ob- 
served, i.e., including when a commensurate second layer 
solid phase is "artificially" stabilized by holding first layer 
atoms at fixed positions, as done in Ref. |20|. Not surpris- 
ingly, no supersolid phase occurs in incommensurate solid 
films, in analogy with what observed in three dimensions. 

In the next section we briefly illustrate our model and 
computational methodology; we then describe our results 
in detail, for the first and second layer. 



II. METHODOLOGY 

Our computer simulations are based on the continuous- 
space worm algorithm^ This methodology has proven 
remarkably effective in large-scale simulations of Bose 
systems, owing to its efficiency in sampling multi-particle 
exchanges, which underlie phenomena such as Bose- 
Einstein Condensation and superfluidity. 

The microscopic model utilized here is standard. 
Specifically, we use the helium Aziz interatomic 
potential ) 22 ' 23 as well as the anisotropic 6-12 graphite- 
helium potential of Carlos and Colef^l used in essentially 
all previous simulation workj 2 ^ Such a potential accounts 
for the corrugation of the graphite substrate, a crucial 
ingredient to describe the variety of registered phases in 
the first layer. Effects of the corrugation of the graphite 
substrate become negligible for successive adlayers; for 
simulations of the second layer we have therefore utilized 
the laterally averaged version of the Carlos-Cole poten- 
tial (see also Ref. [200 . 

Our model is fully three-dimensional, we use standard 
periodic boundary conditions, and simulate systems com- 
prising up to 600 particles, at temperatures as low as 0.2 
K, which is low enough to yield essentially ground state 
estimates. 



III. RESULTS: FIRST LAYER 

In the first layer we confirm the existence of a com- 
mensurate Ci/ 3 solid at coverage p x / 3 = 0.0636A (Fig. 
[2]). At higher coverage, the film enters a domain- wall 
phase (DW), with stripes of the C1/3 solid separated by 
(superheavy) domain walls. 17 

At even higher coverage, we observe a change from 
striped to hexagonal network of (heavy) domain walls. 
The possibility of a transition between these two domain- 
wall types was already raised by Greywalli This network 
becomes denser with increasing coverage, ending with 

a commensurate solid (C7/ 16 ) for Pt/iq — 0.0835A , 
where 7/16 of the adsorption sites are occupied. This 
structure is also found in diffraction experiments of D2 



on graphite^ but had not yet been predicted for Helium. 
Greywall proposed a particular commensurate structure, 

around a coverage 0.820 A , which differs only by 1.8% 
from P7/16- Thus, the signal observed in his heat capacity 
measurements may stem from the C7/16 solid. For den- 
sities above P7/16 we find an incommensurate solid, in 
agreement with experiments and previous calculations. 16 

Next, we determine the layer promotion density p^ , at 
which the second layer starts to become populated. At 
equilibrium, the chemical potentials of the first and the 
second layer are equal, i.e., p\{p\, P2) = ^2(^1,^2), where 
pi and P2 are the densities in the first and second layer 
respectively. Layer promotion occurs at specific value 
of the chemical potential p lp at which the second layer 
density P2 jumps from zero to a finite value p™ m . Thus, 
we can find p 1 / 3 by solving the equation 

^ = m (pi p ) o) = M2(pi p ,pr n ). (1) 

In order to determine p\(pi), we performed simulations 
with a fixed number of particles, and with the first layer 
initialized as a triangular solid. Density scans are ob- 
tained by varying the area of the simulation cell. It is 

de 

Mi(Pi) = e(pi) + pi — , (2) 
dpi 

with e = E/N the energy per particle, where we estimate 
the derivative from a polynomial fit of degree 4 to e(p) as 




FIG. 2: (Color online) Contour plots of the density distribu- 
tion of the 4 Heatoms in the first layer (x-y-plane parallel to 
the substrate) . The black stars mark the minima of the under- 
lying graphite potential (adsorption sites). The distance be- 
tween two adsorption sites is 2.46A. The snapshots are taken 
at different densities in the first layer (see Fig. [IJ: C1/3: com- 
mensurate 1/3 solid, DW: domain-wall phase, Cr/ie- com- 
mensurate 7/16 solid and IC: incommensurate solid. 
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FIG. 3: (Color online) Left plot: Energy per particle versus 
coverage of simulations with a constant number of particles 
(56) and varying volume. The red line is a polynomial fit of 
degree 4. Right plot: The red curve corresponds to the chem- 
ical potential in the first layer, with an error bar given by the 
red dotted lines. The black dashed line indicates the chemical 
potential, at which the second layer becomes populated. The 
error bars are smaller than the symbol sizes. The intersection 
yields the layer promotion density p 1 / 5 = 0.1140 ± 0.0003A 
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FIG. 4: (Color online) Upper panel: 4 He superfluid density. 
Lower panel: peak values of the static structure factor. Data 
shown are for the second layer a temperature T — 0.5 K. 
The system sizes are given as multiples of the first layer unit 
cell. The first layer density is p 1 = 0.1202 A -2 (0.1270 A -2 ) 
in simulations with active (fixed, stars) first layer particles. 
Statistical errors are of the order of symbol sizes. 



shown in Fig. [31 We determine the chemical potential of 
the second layer at the minimum coverage p 2 nm by taking 
the thermalized first layer simulations and switch to a 
grand canonical simulation, where we gradually increase 
(i until the second layer becomes populated. We find ji2 
= —29.6 ± 0.3 K (essentially independent of pi). The 
solution of Eq. |T]) for is given by the intersection of 
the two chemical potential curves in the right panel in 
Fig. [3l Comparing results from different simulations (30 

and 56 particles), we find p^O.1140 ± 0.0003 A" 2 . 

This value is in agreement with previous simula- 
tion results of Whitlock^ who used an effective po- 
tential for second layer particles. It is also compatible 
with neutron diffraction experiments and heat capacity 
measurements; 10 ' 11 ' 13 ' 19 where values in the range 0.112- 
0.115 A - ' were found. 

Grey wall and Busch determined a value of 0.12 A 
from heat capacity measurements^ higher than our re- 
sult. However, in Ref. |l| Greywall pointed out that there 
is an ambiguity in his coverage scale by several percent. 
Thus, our layer promotion density is consistent with their 
value, taking their uncertainty into account. We have 
tested the sensitivity of our result upon deepening the 
attractive well of 4 He-graphite potential by 10%; we find 
Pi = 0.1165 ± 0.0005 A in this case, still lower than 
that of Greywall and Busch. In fact, the value of density 
corresponding to first layer promotion proposed by Grey- 
wall and Busch, is only observed by making the potential 
more attractive by over 20%, a correction which seems 
unlikely, as we discuss below. 



IV. RESULTS: SECOND LAYER 

We now discuss the physics of the second adsorbed he- 
lium layer. Adsorption of successive layers causes a com- 
pression of the first adlayer, experimentally estimated 
between four— and six percent^ In this work, we per- 
formed simulations for the second layer based on differ- 
ent first layer densities, ranging from 0.1164 to 0.1270 

n —2 

A , the latter being the value proposed by Greywall 
and Busch, and assumed by Pierce and Manousakis 
in their PIMC simulations ! 20 ' 27 

On varying the density of the first layer in the above 
range, the physics of the second layer does not change 
qualitatively; typical results are shown in Fig. [4j Specif- 
ically, we only find a superfluid and an incommensurate 
solid phase, separated by a first-order phase transition. 
We do not observe a commensurate solid phase sand- 
wiched between the liquid and the incommensurate crys- 
tal, in contrast to the prediction of Greywall and Busch. 

The 2d equilibrium density of the liquid second layer, 

in the T — > limit, is estimated at 0.046(3) A , essen- 
tially independent of first layer density. In order to give 
an idea of the weakness of the adsorption potential seen 
by second layer atoms, one may note that the above equi- 
librium density is indistinguishable from that of purely 
two-dimensional 4 He j2£ and significantly lower than that 
of a He mono-layer on a lithium substrate (the weakest 
known substrate wetted by 4 He). 30 This liquid film turns 
superfluid at low T, as we established by direct computa- 
tion of the superfluid density pg, based on the well-known 
winding number estimator j21 

9 —2 

Assuming a first layer density of 0.1202 A (cor- 
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FIG. 5: (Color online) Snapshot of the commensurate 7/12 
solid in the second layer, which only appears in simulations 
with fixed first layer particles with a first layer density of 
pi — 0.1270A . The black dots mark the positions of the 
first layer particles (repulsion sites). Second layer particles 
found at the intersections of the thick dotted (yellow) lines 
are located above the middle of three neighboring first layer 
particles. The (yellow) thick line marks a unit cell of the 7/12 
solid. The 4 He superfluid density in this phase is zero. 
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FIG. 6: (Color online) Scaling of the structure factor peak 
with system size A (area) at temperature T — 0.5K. Results 
obtained from active/fixed first layer particles for two differ- 
ent first layer densities pi are shown. The incommensurate 
solid (triangles, IC) is always stable, independent of the first 
layer density. The second layer density is p2 = 0.080 A in 
this example. For all the other data points the coverage cor- 
responds to the one of the commensurate 7/12 solid. Only 
in the case of fixed first layer particles and highest first layer 
density pi = 0.1270A the commensurate solid is stable in 
the thermodynamic limit. 



responding to a 5% compression), the onset value of 
coverage for the occurrence of superfluidity is therefore 
0.166(3) A~ 2 (Fig. [[]). As shown in Fig. [J (crosses and 
triangles), the superfluid density vanishes at a 2d den- 
sity of 0.076 ± 0.002 A , corresponding to a coverage of 
0.196 ± 0.002 A , at which the incommensurate crys- 
tal phase appears. These coverages are altogether con- 
sistent with existing measurements, once experimental 
uncertainties are properly taken into account 4 

In the heat-capacity measurements of Greywall and 
Busch, a peak appears around a coverage of 0.197 A , 
which they associate with a commensurate structure. 
They assumed a compressed first layer density of 0.1270 

A . From the ratio of first and second layer densities 
(w 4/7) , Greywall conjectured a v7 x y7 registered struc- 
ture with one-quarter of the second layer atoms located 
directly above the first layer atoms; this has also been 
proposed for the second layer of 3 He on graphite^ To 
our knowledge, however, no direct experimental evidence 
for the commensurate phase has been reported. One pos- 
sible scenario is that the observed peak in the heat ca- 
pacity is a signature of the incommensurate phase, which 
in our simulations appears at the same coverage. 

In the PIMC simulations of Ref. SO, the 4/7 com- 
mensurate structure was observed (with a different po- 
sitioning of the second layer atoms with respect to the 
first layer), assuming the first layer density proposed by 
Greywall and Busch. However, as mentioned above, in 



these simulations first layer particles were held fixed, i.e. 
their zero-point motion was neglected. On performing 
the same simulation, we also observe a commensurate 
structure (Fig. 2]), albeit at a slightly different filling 
(7/12). It consists of a triangular lattice rotated by an 
angle of 10.89° with respect to the first layer, as shown 
in the snapshot in Fig. [5l The difference in density com- 
pared to the 4/7 filling is only w 2%. The superfluid 
density of such a commensurate phase is zero, i.e., no 
evidence of a possible supersolid phase is found (see Fig. 

As shown in Figs. [?] and O no commensurate phase 
arises if first layer particles are simulated explicitly, even 
if the (relatively high) first layer density of 0.127 A 
utilized in Ref. [2(J is assumed. This is because, due 
to zero-point motion, first layer atoms occupy a larger 
region of space than predicted classically. Thus, second 
layer atoms are slightly pushed away from the first layer 
(see Fig. [7]), which reduces the substrate attraction by 
~ 3 K, enough to de-stabilize the commensurate phase 
(see also Ref. \3^). It should also be noted that, even 
with fixed first layer particles, no commensurate phase is 
observed if a first layer density lower than 0.127 A is 
assumed (Fig. [5J). 

Regardless of whether first layer particles arc held fixed 
or not, at sufficiently high coverage an incommensurate 
crystalline phase forms (see Fig. 0J. Such a phase is also 
not superfluid, consistently with what is now regarded as 
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FIG. 7: (Color online) Comparison between the density pro- 
files of a system with fixed (dashed line) and active (full line) 
first layer particles at a total coverage p = 0.2117 A 2 . In 
the latter case the second layer is displaced by ~ 0.3 A away 
from the substrate due to the zero-point motion of the first 
layer particles. 



a fairly general theoretical statement. 5 

Just like for the issue of first layer promotion, we have 
explored the possibility that a revision of the helium- 
graphite potential might indeed stabilize the commensu- 
rate second layer phase which we do not observe using 
the Carlos- Cole potential as originally proposed in Ref. 
l24l . Simulations with a 10% deeper substrate potential 
also failed to yield a stable commensurate solid, except 

2 —2 

at the very largest first layer density pi = 0.1270 A 
(which is however not compatible with our value for p 1 / 3 ), 
where the 7/12 commensurate solid becomes stable. We 
discuss below whether a quantitative revision of the po- 
tential seems justified, in light of known experimental 
facts. 



V. CONCLUSIONS 

We have carried out a thorough computational study of 
the first and second layers of 4 He adsorbed on a graphite 
substrate, based on the most realistic interaction poten- 
tials currently available, utilized in all previous simula- 
tion work. For the first layer, our results largely confirm 
those of others, and are consistent with experimental re- 
sults. 

Conversely, our simulations show no commensurate 
second layer phase, in contrast to the interpretation of 
heat capacity measurements by Greywall and Buschii^. 
Rather, we only observed a liquid layer (superfluid at low 
temperature), which crystallizes at high density to form 
an incommensurate solid. The physical behavior of the 
second layer is thus very close to that of a purely two- 
dimensional system, also looking at the narrow spread 



of the 4 He density in the perpendicular direction and at 
the non-existent overlap of first and second layer parti- 
cles shown in Fig. [7j pointing to the absence of inter- 
layer quantum exchanges. It is also worth pointing out 
the similarity between the physics of the second layer of 
4 He on graphite and that of a mono-layer adsorbed on a 
lithium substratej 30 i 34 

Our prediction of no commensurate phase, is at 
variance with previous numerical work by Pierce and 
Manousakis, who observed it instead in their PIMC sim- 
ulations. It need be stressed, however, that there is no 
significant numerical disagreement between our results 
and those of Ref. [2(| Our different conclusion directly 
stems from the fact that, unlike Pierce and Manousakis, 
we did not keep first layer particles fixed, but rather sim- 
ulated their zero-point motion explicitly. This fact alone 
accounts for (most of) the difference between the physi- 
cal outcomes of the two studies. It is worth noting that 
the possible importance of the role of quantum fluctua- 
tions of inner layer particles had already been suggested 
by other authors^ 3 - 

The absence of a commensurate phase suggests that 
either the experimental data have to be reinterpreted, or 
the microscopic model adopted so far, chiefly the helium- 
graphite potential, might have to undergo significant re- 
visions. We have studied this scenario to same length, 
in this work, notably by rendering the attractive well of 
the potential deeper. If the attraction is increased by 
some 10% (while fully retaining in the simulation the 
zero-point motion of first layer 4 He atoms) a commen- 
surate crystalline phase is observed only if a density of 
0.127 A is assumed for the first layer, i.e., the value 
proposed by Greywall and Busch. This value is consis- 
tent with the first layer promotion density computed in 
this work with such a revised potential (i.e., 10% more 
attractive) , only on assuming a compression of the inner 
adlayer of some 9% upon adsorption of successive layers. 
This degree of compression is substantially above that 
estimated by most experimental studies. 

In any case, a revision of such a quantitative degree 
of the helium-graphite potential proposed by Carlos and 
Cole seems doubtful, in view of the good agreement with 
experimental data presented in Ref. l24l 

Finally, no finite superfluid signal is seen in any of 
the crystalline phases observed, either the incommensu- 
rate or the commensurate (the latter, as explained above, 
merely occurs as the result of treating inner layer parti- 
cles as fixed in the simulation). Thus, we conclude that 
this system is no realistic candidate for the observation 
of a supersolid phase. 
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